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ABSTRACT 

The problem of interaction of the rotating magnetic field, frozen to a star, with a thin well 
conducting accretion disk is solved exactly. It is shown that a disk pushes the magnetic field 
lines towards a star, compressing the stellar dipole magnetic field. At the point of corotation, 
where the Keplerian rotation frequency coincides with the frequency of the stellar rotation, 
the loop of the electric current appears. The electric currents flow in the magnetosphere only 
along two particular magnetic surfaces, which connect the corotation region and the inner 
C"| . edge of a disk with the stellar surface. It is shown that the closed current surface encloses the 

Q_f magnetosphere. Rotation of a disk is stopped at some distance from the stellar surface, which 

is 0.55 of the corotation radius. Accretion from a disk spins up the stellar rotation. The angular 
momentum transferred to the star is determined. 
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1 INTRODUCTION 

The problem of interaction of an accretion disk with a magne- 
tized star is very important for understanding of compact X-ray 
sources. The energy release for a variable X-ray sources and X- 
ray pulsars is due to the accretion of the matter from a companion 
star onto the surface of a compact star, namely, on a neutron star. 
Falling onto a neutron star, a particle gains the energy per unit mass 
equal to the value of the gravitational potential on the stellar sur- 
face ip = 0.15c 2 (M s /M©)(.R ii /10fcm)~ 1 . where M s and R B are 
stellar mass and radius, respectively, and Mq is the solar mass. 
The energy release depends on how the accretion is realized. The 
magnetic field, which is frozen to a neutron star, is high enough to 
prevent a free falling of the matter on a star. The naive point of view 
is that the accretion is free outside the so called Alfven surface, and 
then continues along the magnetic field lines down to the polar caps 
of a star (Lamb 1989). In a standard approach Alfven surface is the 
surface where the magnetic field pressure B 2 /Sty is equal to the 
pressure of the accreting matter. But here it is not clear what is the 
pressure of the accreting gas. Is it the total dynamic pressure gv /2 
or its radial part only, gv 2 /2 ? The Alfven radius is estimated for 
the spherical accretion, v = v r , but then is applied to the disk ac- 
cretion. But in this case the gas velocity v comes mainly from the 
velocity of disk rotation and is of the order of Keplerian veloc- 
ity vk , which is much larger than the radial velocity. The gradient 
of the magnetic pressure in the azimuthal direction is small and 
therefore does not oppose significantly the disk rotation. Besides it 
is clear that the disk must compress the magnetic field within the 
disk inner edge. These features do not help to answer the question: 
where is the real position of the Alfven surface. 

An accreting gas is a plasma consisting of ions and electrons, 
which have a quite different Larmor radii in a magnetic field. Mov- 
ing across the magnetic field, the plasma gets polarized and creates 



the large scale electric field E. If this field is E = v x B/c 
then the plasma continues to drift in the crossed B,E fields at 
the same velocity v not slowing down. The real value of E de- 
pends on the boundary conditions on a disk and a star surface, 
on electric currents in a plasma and on conditions of their clos- 
ing. The important characteristics are the electric conductivity of 
a disk plasma and the neutron star conductivity in the surface lay- 
ers. The ionized plasma of an accretion disk has the conductivity 
a = 10 13 (T e /leV) 3/2 (A/10)~ 1 s~\ which is high enough to 
consider a disk as an ideal conductor. Here T e is the temperature of 
electrons in a disk, which is higher than 10eV; A is the Coulomb 
logarithm (A ~ 20). At such a high conductivity a the width of 
the skin layer X s k = (rc 2 /<r) 1/ ' 2 is less than the disk width H. 
The value of r is the characteristic time of turbulent motion in the 
Q-disk, r = H/vk- The condition a ^> c 2 /Hvk is well fulfilled 
in the inner parts of a disk. As for the conductivity of surface layers 
of a neutron star ctats, then it is as high as 10 21 s _1 (Blandford et 
al. 1983; Potekhin 1999). Therefore, ctats 2> <r, and a neutron star 
can be considered as ideal conductor too. 

As ideal conductor, a disk tends to exclude the stellar magnetic 
field, pushing it toward a star. Heavy disk ions tend to rotate with 
the Keplerian velocity, but magnetized electrons are frozen to the 
magnetic field lines, which rotate with the angular velocity of a star 
lj s . Thus, the point p — p c where a disk corotates together with a 
star, uj 3 p c = (GM s /p c ) 1/2 , 
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is the point beginning from which (p < p c ) motions of ions and 
electrons differ essentially. Here G is the gravitational constant and 
P 3 — 2n /us is the period of the stellar rotation. The region of coro- 
tation is namely the region where the interaction of an accreting 
disk with the stellar magnetic field begins. The point of corotation 
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is inside the stellar light cylinder radius Rl = c/oj s for all rotating 
neutron stars: P s > 3 • 1CT 5 (M s /M Q )s. On the other hand, the 
corotation radius is larger than the neutron star radius R s . Because 
of that we will consider the undisturbed magnetic field of a neutron 
star as dipolar. Also we assume that the axis of a dipole is parallel 
to the neutron star rotation axis. 

The paper is organized as follows. In section 2 we describe 
an arbitrary axisymmetric magnetic field. In section 3 we find the 
motion of ions in a disk. Section 4 describes electric currents in 
the stellar magnetosphere. In section 5 we find the structure of the 
magnetic field. Finally, in Section 6 we find the torque acting on a 
star and discuss the obtained results. 



2 AXISYMMETRIC MAGNETIC FIELD 

Since our aim is to find the magnetic field of the magnetosphere of 
a star, we initially introduce convenient variables, which describe it 
more simply. For an axisymmetric magnetosphere, when the angle 
between the axis of the magnetic dipole and axis of the rotation of 
a star equals zero, it is useful to introduce the flux of the poloidal 
magnetic field f(p, z). Here p and z are the cylindrical coordinates. 
Then the components of the magnetic field are 
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The function g(p, z) describes the toroidal magnetic field produced 
by the poloidal electric currents flowing in a stellar magnetosphere. 
The relation / = constant is the equation for the magnetic surface 
on which the magnetic field lines are lying. It is also convenient to 
describe magnetic field lines of the poloidal magnetic field not only 
as / = /(Pi 2 )> but also p — p(z, /). The poloidal magnetic field 
is 
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For the dipole magnetic field the magnetic flux f(p, z) is 

B a R\p 2 

U (p 2 + 2 2 ) 3 / 2 ' 

where B s is the magnitude of the magnetic field on the surface of 
a star on its equator, B s = \B z (r = R s )\ eq . For our problem it 
is important to know the value of magnetic flux in the equatorial 
plane, 2 = 0, fo(p) = f(p,z = 0). Then the vertical magnetic 
field on the equator is B z (z = 0) = p~ 1 dfo/dp. 

It is useful to measure the magnetic field in the units of B s 
and the distances p, z in terms of star radius R s ; then the units for 
the magnetic flux will be B S R 2 S . The relations (1) are the same in 
these dimensionless variables. For a dipole field, the dimensionless 
flux is f d = p' 2 /(p' 2 + z' 2 f'\ f' 0d = l/p>. 

In the star magnetosphere, electric currents flow along the 
magnetic field lines, j = aB, where a(r) is an arbitrary scalar. 
In this case, the poloidal components of the Maxwell equation, 
V x B = 47rj/c, result in the relations 

4tt d/' 

These relations mean that the toroidal magnetic field and the elec- 
tric currents in an axisymmetric magnetosphere are functions of 
the poloidal magnetic flux, /. The important case is the absence of 
toroidal volume electric currents, — 0. Then (V x B)^ = 0, 
and the equation for the magnetic flux / is the Laplace type 
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The function / is neither even nor odd with respect to the coordi- 
nate z due to the electric currents flowing in accretion disk in the 
equatorial plane. Thus, we expand the function / over the exponen- 
tial functions cxp( — Xz), for z > 0. The solution of the equation 
(2) is 



f(p,z)= f 
Jo 



exp(—Xz)XpJ 1 (Xp)<p(X)dX; z > 0, 
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where Ji is Bessel function of the first kind. An arbitrary func- 
tion tp(X) in this integral is defined by boundary conditions for the 
equation (2). For the dipole magnetic field function ip is ip(X) = 1 
(in dimensionless units). The deviation of ip from the unity de- 
scribes the distortion of the stellar magnetospheric magnetic field 
from dipolar. It is useful to express the function ip through the value 
of the magnetic flux on the boundary z = 0, /o(p). Because of 
fo(p) = J °° XpJi(Xp)cf)(X)dX, using the inverse transformation, 
we get 
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The equations (3) and (4) determine the magnetic field through the 
boundary condition /o(p). Of course, we also take into account an- 
other boundaries: infinity and the axis p = 0. The quantity f(p, z) 
is finite there except the point p 2 + z 2 — > 0, where / — > f d . 



3 THE ACCRETION DISK 

The disk is produced from the matter falling from a companion star 
to a neutron star. The rate of accretion is of M. The centrifugal 
force prevents a disk matter from spreading over a magnetosphere, 
and this gives rise to an equatorial plasma disk. The width of this 
disk is small compared with its radial size. The strong magnetic 
field of a neutron star forces the disk plasma to rotate with the 
angular velocity of a star. We will describe motions of ions and 
electrons separately because their motions are quite different in the 
strong magnetic field of a star. Disk ions have radial, v p , and az- 
imuthal velocities, v$. We consider a stationary and axisymmetric 
star magnetosphere. Due to axial symmetry, all quantities are inde- 
pendent of the azimuthal angle (f>. Then, the continuity equation for 



\h {pniVp) = ^i 5{p ~ pout)s{zl 



(5) 



where ni is ion density, p out is the outer radius of the disk, p out ^> 
p c . Integrating Equation (5) over p and z, we obtain 
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Here S is the surface disk density, E = J nidz. Let us note the 
very important fact that the divergence of the ion flux at p < p out 
equals zero. 

Now we write the equation of motion of ions 



m;(v • V)v = QiE + — v x B + F. 



(7) 



The star acts on ions by gravitational force, F p = —miGM s /p 2 , 
which, neglecting other forces, must result in the rotation of ions 
with the Keplerian velocity, vk = {GM S / p) 1 ^ 2 . The p-component 
of Equation (7) gives 
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We can see from Equation (8) that in the absence of the electro- 
magnetic fields B z and E p , the ions rotate with the Keplerian ve- 
locity, v<f, — vk- In contrast, in a strong magnetic field B z , ions 
have the electric drift velocity, — —cE p /B z . For corotation, i.e. 
the rigid rotation with the angular velocity of the star, v$ = uj s p, 
the plasma in the disk must be polarized to create the radial elec- 
tric field, E p = —u s pB z /c. The corotation velocity is less than 
the Keplerian velocity at the distance p < p c . Because the Keple- 
rian velocity goes down with increasing distance p, the corotation 
of a plasma disk means that the centrifugal force in the disk be- 
comes less than the gravitation force at p < p c , and ions must 
move inwards in the radial direction. But this is not easy, due to 
the conservation of angular momentum. That is described by the 
(^-component of the equation (7), 
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Let us note that there is no the toroidal electric field, E^ — 0, 
because under stationary conditions the electric field must be po- 
tential. Because the ion radial velocity, v p , is not equal to zero (see 
Equation (6)), the solution of Equation (9) is 



pvs H —fo(p) = const. 
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Equation (10) is the conservation law of the total angular mo- 
mentum of an ion in magnetic field. But the dimesionless co- 
efficient before the second term in Equation (10), written as 
(qiB /cmiOj s )(uj 3 f(,l B), is much greater than unity. It is propor- 
tional to the ratio of the ion cyclotron frequency of rotation in the 
magnetic field of a star, uj ci = qiB s /cnii, to the frequency of ro- 
tation of star, aj s . We denote this ratio Q, c = uj ci / ui s . This is the 
cyclotron frequency of ions in natural units of the stellar rotation 
frequency 
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For the dipole magnetic field, the terms in the left hand side of 
Equation (10) become of the same order only at large distances 
pjR s — n^ 3 . Angular momentum conservation prohibits ion ra- 
dial motion in such a strong magnetic field as that of a star if 
we assume it is (close to) dipolar one. The only possibility to al- 
low ions to move radially is an expulsion of magnetic field from 
the disk. Only if the magnetic flux changes slowly in the disk, 
A/o ~ B s p fic , does the motion of ions in radial direction be- 
come possible. But a small variation of the magnetic flux fo with 
the radial distance p means a small value of the vertical magnetic 
field B z in the equatorial plane (see formula (1)). 

As for electrons, the parameter S7 C for them is at least 10 3 
times greater than for ions. The plasma electrons are strongly mag- 
netized and can move only along a magnetic surface / = const. 
They generate a current stream at / = f c = fo(p = p c ) towards 
a star. Thus, the disk plasma is separated at p < p c : ions create the 
plasma disk in the equatorial plane, but electrons create the electric 
current J c on the magnetic surface f = f c . This does not mean that 
the plasma disk is positively charged, it is neutral. The electrons 
from the 'sea' of electrons of a stellar magnetosphere neutralize 
any electric charge. 



4 MAGNETOSPHERIC AND STAR SURFACE 
CURRENTS 

We consider the magnetospheric plasma which has an infinite con- 
ductivity. Under stationary conditions it can rotate around the z axis 
with the angular velocity ui(p, z). In the frame rotating with this ve- 
locity the electric field E' is equal to zero. So, the electric field in 
the laboratory frame is 
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c dp 



The electric field must be potential E — — V*, where $ is the 
electric potential. This means that the frequency of rotation w(p, z) 
is only function of the magnetic flux f,u = w(/), and the electric 
potential is * = J f u(f')df'/c Because the conductivity of a star 
surface is very high, w(/) = w s . 

The magnetospheric current J c flows from a star surface to 
a disk along the magnetic surface / = f c . This current is closed 
by the star surface current J s . The surface current, J s , flowing on 
the stellar surface is connected with the volume magnetospheric 
current, j, by the continuity equation 

V s ■ Js = jn • 

The divergence in this equation is taken along the surface where 
surface currents flow, and j„ is the normal (to this surface) compo- 
nent of the volume current j. The analogous equation connects the 
ion surface current of the disk J d with the magnetospheric volume 
current j. But the divergence of the ion surface current, as we have 
already emphasized (see equation (6)), equals zero. The electrons 
of the disk can only rotate with the angular velocity of the magne- 
tosphere u(f), and the divergence of their flux is also equal zero. 
Thus, the volume magnetospheric current j is absent except on the 
particular magnetic surfaces / = f c and / = f s . The value of 
/ = f s is the magnetic flux where the disk stops rotation, f B > f c . 
The total current / produced by the accretion flow is 
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One half of it flows in the north hemisphere and another half flows 
in the south hemisphere. The magnetospheric currents flow along 
the magnetic field on two magnetic surfaces 
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These currents are closed by the surface current flowing in the stel- 
lar polar region. 

The magnetospheric electric currents j s , j c produce in the re- 
gion f 3 > f > f c the toroidal magnetic field B^, 
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Then the (^-component of the Maxwell equation V x B = 47rj / c 
gives the equation for the poloidal magnetic flux 
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+ ^ = \W-fs)-5{f-f c )\. (14) 



The right hand side of the equation (14) is not equal to zero only 
on the magnetic surfaces / = f s ,f c , where the magnetospheric 
electric currents flow along the magnetic field. Due to the toroidal 
component of the magnetic field the longitudinal currents have 
a toroidal component, which distorts the poloidal magnetic field. 
Outside the surfaces f s , f c , Equation (14) coincides with Equation 
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(2) and its solution is given by relations (3,4). However, the fields 
on both sides of the discontinuities are different. The connection 
between the magnetic fields at / > f s and / < f s , and also be- 
tween / > f c and / < f c , can be obtained by integrating the 
equation (14) over / in the small region near surfaces / = f s and 
f ~ fc respectively. For that we need to go over to another inde- 
pendent variables, (z, /), and consider the magnetic surfaces to be 
determined by the relation p = p(z, /). Equation (14) becomes 
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We take into account that on both sides of discontinuity the shapes 
of the magnetic surfaces are the same 
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Then integrating the equation (15) we obtain 
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(16) 



This is the natural condition of equality of the magnetic pressures 
on both sides of a discontinuity, which follows exactly from our 
equations. 




5 STRUCTURE OF MAGNETIC FIELD 

Here and below we will use the dimensionless variables 

J = cj/cj s , r' = r/R s , v' = v/R s oj s , B' = B/B s , 
f = f/B s R 2 s , £' = 2-KrrnUsRlT.IM, J' = 4nR s J/I. (17) 

Further on we will omit the index /. There appear several dimen- 
sionless parameters. The ratio of the ion cyclotron frequency to the 
frequency of star rotation, mentioned alrready in the section 3, is 
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The second parameter is the magnetic field, created by the current 
7/2(11) flowing on the star surface, in the units of B s 



B : 



qiM 



cmtRsBs 



B = 60 



(19) 



M 



10- 10 M e /y J U0 12 G 



B s 



R, 



10 6 cm 



Moving toward a star a disk pulls out the star magnetic field 
lines. Vertical component of the magnetic field B z decreases and 
the ions can move in the radial direction. But the electrons remain 
magnetized. They must move along magnetic field surfaces only. 
The single possibility for them to change their motion from the 
Keplerian rotation to the corotation with a star is to do that at the 
point of corotation p — p c . Thus, a disk banks all dipole magnetic 
field lines crossing the equatorial plane at p > p c into the region 
p < p c - This must had happened because magnetic field lines were 
froze initially to the ideal conducting matter of a disk and had to 
move together with it. We conclue therefore that the stellar magne- 
tosphere under the action of an ideal accretion disk is confined to 
the region p < p c l Thus, the magnetic field flux at the point p = p c 
equals zero, f c = 0. The magnetic field line / = just comes from 



Figure 1. Topology of the magnetic field and currents. Thick lines are elec- 
tric currents. 



the magnetic pole of a star. The topology of magnetic field line is 
drawn on the figure 1 . 

The velocity of rotation of plasma disk ions is given by Eq. 

(10) 



_ pc_ _ n e /o(p) 
p p 



(20) 



Here we consider that at the point p = p c the velocity of ion 
rotation is equal to the Keplerian velocity vk- Eq. (20) fixes the 
value of magnetic flux at the stopping of the disk rotation — 0, 
f s — p^/Qc- The electrons of the disk rotate with the angular ve- 
locity of rotation of the magnetic field lines uj(f) = oj s , however. 
Thus, the plasma disk produces the surface toroidal electric current 



J 4, = 2SE(^ - p), 



(21) 



which gives the discontinuity of the radial magnetic field on the 
surface z = 0. 



B P \ z=z+0 = - B p \ z=z _ = BE(«* - p). 



(22) 



The value of (v$ — p) is a function of p, fo(p) and is determined 
by the relations (20) 
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It is seen from (21,23) that the value of the radial magnetic field up 
to the stopping point p = p s is negative, B p (p, z — +0) < 0. At 
the point p = p c the radial magnetic field in the equatorial plane 
is equal zero, and the magnetic field line here is perpendicular to a 
disk. Then magnetic field lines bend toward a star. 

The value of the radial magnetic field B p (p,z — z + 0) de- 
termines the derivative of the magnetic flux over z on the equator 
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Let us find this derivative as a function of fo (p). From the equations 
(3,4) we obtain 
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The second integral is known (Oberhettinger 1972), and we get 
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Here K(x) is the complete elliptic integral of the first kind. It 
equals n/2 at x — 0, but has the logarithmic singularity at x — > 1, 
K(x) — > ln[16/(l — x 2 )]/2. Equation (25) must be supplemented 
with the condition of absence of the radial magnetic field in the 
region p < p s , dfo/dz — 0. Let us note that substitution of the 
dipole magnetic field, fo = 1/p, into Eq. (25) satisfies this con- 
dition for all p. Also the square of the vertical magnetic field B 2 
behind the stopping point, p — p s — 0, must be equal to the sum 
of the square of the the radial one B 2 (B z « B p ) and the square 
of the toroidal one B\ — B 2 / p 2 before the stopping point (see 



the condition (16)). As a result we get B z 



-B(l + P i/vi(p 



Ps)) 1 ^ 2 1 ps- We see that the magnetic field is strongly compressed 
here by a disk, \B Z \ >> [Bdip| = Ps~ 3 - The equations (24-26) 
close the system of equations which determine the function /o(p). 
We add also the energy conservation law following from the equa- 
tions (8,9) 
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Qcfo(p) = const(p). 



(26) 



The particle energy decreases due to the work of the radial elec- 
tric field in the disk on ions. This electric field is generated in the 
disk by the rotating magnetic field. Considering that the radial disk 
velocity v p is smaller than the rotation velocity v$ — vk at the 
point p — p c , we find the radial velocity of a disk at any point 
ps < p < p c 

„3 
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pc - v<t> - 2n e / (p). 



The radial velocity v p defines the surface density E(p) — 
l/(p|« p |). The radial velocity v p at the stopping point p = p s does 
not equal zero and is negative 



v P (p = ps) = -pc 
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At the point p — p s disk ions begin to move along the magnetic 
surface / = f s toward the star surface to its polar cap. Particles 
without the mechanical angular momentum (v^ — 0) fall down to 
a star freely but along curved fixed magnetic lines. They get the 
kinetic energy v 2 /2 = p\ — p 2 on the stellar surface which corre- 
sponds to fall down from the height of the corotation point p c . 

The function fa(p) on all axis < p < oo is different in three 



regions: f = f^'(p) at p < p 3 ; f = f^'(p) at p s < p < p, 
fo = at p > p c . Thus, the equation (25) becomes 



(2), 



i:l::^.e , (PMdx}=o,p< P s 



(27) 
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Figure 2. Dimentionless poloidal magnetic flux in the equatorial plane, fo, 
versus distance p from the z axis. Four regimes of the p-dependance are 
distinguished: A R s < p < p B — Ap (Ap << p s ) where the dipolar for- 
mula is valid, /o(p) = 1/p; B Narrow region p s — Ap < p < p s where 
the poloidal field lines are strongly compressed by the disk, and /o (p) de- 
creases rapidly from l/p s down to/ s <Cl/p s ; C p s < p < pc - Here, 
/o(p) is small and decreases from f s to /o(Pc) = 0. D For p > p c , 
fo(p) = 0. For details see section 5. 



Ps < p < pc- (28) 



B 



(2) 

The value of the magnetic flux /q is small. At the point p = p s 

(2) 

it equals /q = f s « l/p s , then it decreases to zero at p = p c . 
On the contrary, the function /^'(p) increases strongly from the 
value /q 1 ' = f s at p — p s to the dipole one l/p s at p < p s due 
to the large derivative df^ /dp = -6(1 + p 2 /v 2 p ) 1/2 . The width 
of this transition region Ap is small, Ap ~ l/(p s B) << p s - Here 
the magnetic field is large, it is strongly compressed by the disk. 
Then, approaching to a star, the magnetic field becomes dipolar 
/q 1 ' = 1/p. The function /o(p) is plotted in the figure 2. 

The position of the stopping point p s can be determined from 
the equation (28). At the corotation point p c the right hand side in 
(28) equals zero. The main contribution to the l.h.s. comes from 
the dipole part of fo, /q 1 ' = 1/p. Substituting this dependence, 
integrating, differentiating and equating the result to zero at point 
Pc we obtain the equation 



E\ P - 
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K\ P - 



Pc J \ Pc 



Here E(x) is the complete elliptic function of the second kind. The 
root of this equation, except the trivial one, p s = 0, is p s = 0.55p c . 
The point p = p s can be called the Alfven point. But it corresponds 
to the stopping of the disk rotation, not of the radial motion. And 
its position is fixed by the corotation radius and does not depend on 
the star magnetic field B 3 . That is due to the strong compression of 
the dipole magnetic field by a disk in the region p < p s . 
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6 DISCUSSION 

We have shown that the structure of the rotating neutron star mag- 
netic field, interacting with a thin accretion disk, can be exactly 
solved, and we obtained an analytical solution for this field. To 
achieve this we have used ideal approximations: an axisymmet- 
ric, well conducting magnetosphere with a thin plasma disk. These 
conditions are adequate to the reality. We have demonstrated that 
a disk compresses the stellar magnetic field, pushing it towards a 
star. The magnetic field is compressed in the region p < p s , where 
p s is the inner edge of the disk where it stops rotation in the lab- 
oratory frame, p s = 0.55p c . After that the matter of a disk, with 
zero mechanical angular momentum, falls down onto a star along 
the magnetic field surface / = / s . The loop of the electric current 
starts near the point of corotation p — p c . The electric current / 
(11) flows along the disk at p s < p < p c , then half of it, 1/2, 
flows along the magnetic surface / = f s to the stellar surface , 
then along the stellar surface in the polar region, and finally re- 
turns to a disk along the magnetic surface / = f c = 0. There 
is a symmetric current structure in the second hemisphere. Let us 
emphasize that the magnetosphere exists only inside the magnetic 
surface / = f c — 0. It is maintained by the pinching Ampere force 
fa = je x / c acting on the electron current j e flowing outside a 
star along the magnetic surface f c . This force just balances the gra- 
dient of the magnetic pressure X7B 2 /8tt inside the magnetosphere 
boundary. 

The potential difference * created by the rotating magnetic 
field is (Landau et al. 1984) 



The electric current J, flowing along the star surface perpendicular 
the magnetic field, creates the Ampere force and spins up the star 
rotation. The torque K is 

K = Jtf/w a = Mp 2 c Lo s . (29) 
It produces the acceleration of the star rotation P s 

-^•10- 13 (Tf) 7/3 (To-4^) 

(ma 2/3 ( u r 1 
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Here I s is the moment inertia of a star. This value of P s corresponds 
to the observed one for the close X-ray binaries (Rappaport & Joss 
1983). We see that the total angular momentum, which has a disk at 
the point of corotation, is transferred to a star. The answer does not 
depend on the strength of the stellar magnetic field B s . It is due to 
the strong magnetization of a disk plasma. The magnetization pa- 
rameter fi c (18) is much greater than unity. And no matter whether 
it is very large (10 12 ), or not so much (10 5 ). The role of the stellar 
magnetic field consists in stopping of the disk rotation, i.e. trans- 
formation of the mechanical angular momentum of particles to the 
electromagnetic momentum. The angular momentum of a charged 
particle in a magnetic field consists of two parts: the mechanical 
one and the electromagnetic one (see the formula (10)). And ap- 
proaching a star, a particle loses its mechanical part and increases 
the electromagnetic one. Therefore, if a matter element falls down 
onto a star, M > 0, then it transfers its angular momentum to a star, 
K > 0. No matter whether it is the mechanical angular momentum 
or the electromagnetic one. Thus, under the accretion, there must 
not be the negative part in the torque, which can be attributed to the 
electromagnetic torque (Ghosh et al. 1977). Appearance of the stel- 
lar spin down is connected with assumptions that, first, the toroidal 



magnetic field changes its direction at some value of the radius, 
and, second, the value of the poloidal field B z in the disk does 
not differ strongly from the dipole one (for example, see Kluzniak 
& Rappaport 2007). The toroidal magnetic field is determined by 
poloidal currents, and the zero value of it at some distance means 
the closure of currents. At larger distance the field B^ is equal to 
zero. Indeed, we see that the toroidal field exists only in the region 
p s < p < p c and B<f, > (13). 

Let us note that at the small value of accretion rate M it is 
impossible to construct the self consistent solution for the stellar 
magnetosphere. When the magnetic field at the point p s becomes 
less than the dipole one at this point it means that there appears 
the decompression of the magnetic field. This is impossible for the 
plasma moving toward a star because the magnetic field is frozen to 
the well conducting matter. Formally, in this case it is impossible 
to match the solution for the magnetic flux fo(p) with the dipole 
field. This takes the place when \B z (p = p s )\ < p~ 3 . Because 
\B Z \ = B(l + p 2 s /v 2 p ) 1/2 /p a ,p a = 0.55p c andv p = -0.8p c , this 
condition is B < 2p~ 2 - It means that the expression (29) for the 
torque is valid when M > M m in , 



qi Pc 

Here p s is the stellar magnetic moment, p s = B s Rl- For M < 
Mmin a disk does not penetrate into the star magnetosphere. 

In some sense, the stellar magnetic field does not prevent the 
accretion of a ionized matter onto a star, but, on the contrary, helps 
it. Without magnetic field an ideal matter element can not accrete 
due to the angular momentum conservation. We need to introduce 
an anomalous viscosity represented by a parameter in the a-disks. 
In a magnetic field a charged particle possesses two parts of the an- 
gular momentum, the mechanical one and the electromagnetic one. 
In the language of magneto hydrodynamics, we have the viscous 
stress tensor and the magnetic stress tensor, which is proportional 
to the product B^B Z . Appearance of the toroidal magnetic field is 
impossible without poloidal electric currents, they must be closed 
on the stellar surface, transforming the mechanical angular momen- 
tum of a disk into the stellar angular momentum. 

It is clear that if the inner edge of a disk p d is placed further 
than the the corotation point, then it will extract the angular mo- 
mentum from a star. In this case the centrifugal force is larger than 
the gravitation force, and a disk must be pushed outward from a 
star. This is the propeller regime. It is stationary, or quasi station- 
ary process, if we have a source of the matter at the inner edge of a 
disk. Such situation takes place in the Jupiter magnetosphere where 
the Jupiter's satellite Io supplies a gas to a plasma disk at the point 
pd > Pc- The electric current in the current loop has the direction 
opposite to that in the case of accretion (Istomin 2005). If there is 
no gas source in the stellar magnetosphere then a disk is ejected out 
to the light cylinder taking the angular momentum A J from a star 
of the order of 

AJ~ M d R 2 L u} 3 . 

Here Ma is the mass of a disk and Rl is the radius of the light 
cylinder, Rl = c/cu s . 

Up to now we considered only the situation when a disk and 
a star rotate in the same direction (corotating disk). There exists 
the interesting case when an accretion disk rotates in the oppo- 
site direction than a star. A corotating disk had the same sign of 
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the angular momentum as a star. The disk matter falling onto star 
transmits its angular momentum to a star. Then, the torque is pos- 
itive. If the disk rotates in the opposite direction than a star (coun- 
terrotating disk), then, matter falling onto the stellar surface, will 
spin down the stellar rotation. The torque becomes negative. For 
the countercotating disk it is impossible to build up the stationary 
self consistent solution like we did for the coratation case. A star 
must destroy the inner parts of a disk approaching too close to it. 
It tends to re-rotate the inner part of a disk to the direction of its 
own rotation. In the other words, the magnetosphere-disk coupling 
tends to enforce corotation to the inner layers of the disk. The ac- 
cretion becomes more quasispherical than disk-like. We can esti- 
mate with a good accuracy the torque acting on a star in this case. 
And the exact solution obtained above will help us. We saw that 
to reduce to zero its mechanical angular momentum, a falling par- 
ticle must deflect from its position on the magnetic field surface 
/ by the value of A/ = p m v m /£lc (see formula (20)). Here val- 
ues p m and v m are the radial position and the azimuthal velocity 
of a particle on the magnetospheric boundary, respectively. For the 
rotating disk, the boundary of the magnetosphere is the corotation 
radius, p — p c . For the counterrotation there is no a corotation 
point. The magnetospheric boundary for the quasi spherical accre- 
tion must be close to the Alfven radius, p m ~ pA- Falling onto 
a star, ions move almost along the magnetic field lines deflecting 
by a small value of A/. Electrons will deflect also, but their mag- 
netization parameter Q. c is much larger than that for ions, and we 
can neglect electron deflection. Thus, ions and electrons will be on 
the stellar surface at different magnetic surfaces separated by the 
value of A/. The current 1(11) will flow along the surface, it will 
produce the torque K (29). For the counterrotation of a star and a 
disk the value of A/ will be negative, and the electric current on 
the stellar surface will have the opposite direction than for the case 
of corotation. As a result the torque acting on a star will be negative 
and equal to K — —Mv m p m . That is the exact relation. But the 
position of the magnetosphere boundary is not known exactly. If 
we do not take into account the magnetic field compression on the 
boundary it is given by the traditional formula for the Alfven ra- 
dius p m = pA = (p'i/2 1 / 2 M pc^Us) 2 / 7 (see, for example, Arons 
1993). Because for the accretion onto a star the radius of the mag- 
netosphere p m must be smaller than the corotation radius p c , the 
deceleration torque will be always less than the acceleration torque 
at the same value of the accretion rate M. Thus, the observed stel- 
lar spin down with noticeable M can be a result of the accretion 
of the matter from a counter rotating disk. This is quite possible 
because the disk formation depends on a companion star, on the 
stellar orbital motion and many other reasons. 
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In conclusion we compare the results obtained here with the 
numerical calculations. The numerical simulations were done in the 
frame of the MHD approximation by Romanova et al. 201 1, 2012. 
Though our treatment have been done in the frame of two-fluid 
hydrodynamics, the similar features were seen in numerical single- 
fluid simulations: the compression of the stellar dipole magnetic 
field by a disk, the falling of the disk matter along fixed magnetic 
field lines, the pushing out of magnetic field lines from a disk, and 
the proportionality of the torque K, acting on a star, to the accretion 
mass rate M. 



